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The nucleon form factors of the energy-momentum tensor are studied in nuclear medium in the 
framework of the in-medium modified Skyrme model. We obtain a negative D-term, in agreement 
with results from other approaches, and find that medium effects make the value of d\ more negative. 
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1. Understanding the structure of the nucleon has been one of fundamental issues well over decades. While the 
scalar, vector and axial-vector properties of the nucleon have been studied extensively and comprehended to a great 
extent, its energy- momentum tensor (EMT) form factors have been drawn to attention only quite recently, long after 
Pagcls proposed them pj. The reason can be found in the fact that it is very difficult to get direct access to these 
form factors experimentally. However, the generalized parton distributions (GPDs) which are accessible via 

hard exclusive reactions [6H15| , make it possible to extract information on the EMT form factors of the nucleon. In 
particular, certain Mellin moments of the GPDs can be expressed in terms of the EMT form factors [H, [l6|, E3- 

The nucleon matrix element of the total symmetric EMT are parameterized by three form factors as follows|16l. [ItJ 



(pi\f^(0)\ P )=u(pi, s') 



M 2 (t) ^ + J{t) \ + dl(t) 



M N y ' 2M N v ' 5M N 



u(p, s) , (1) 



where P = (p + p')/2, A = (p' — p) and t = A 2 . Mjv is the nucleon mass, and u(p, s) denotes the nucleon spinor 
with the polarization vector s defined such that it is given as (0, s) in the rest frame in which s designates the axis of 
the spin quantization. The recent interest about the EMT form factors was stimulated by the fact that it is possible 
to define in QCD and to access in experiment the separate quark and gluon contributions to the form factors. The 
only presently known non-trivial piece of information is the decomposition of M^it) at the zero-momentum transfer, 
which reveals that about 1/2 of the momentum of a fast moving nucleon is carried by quarks, and the other half by 
gluons. J(t) provides analogous information on how the total angular momentum of quarks and gluons makes up the 
nucleon spin, but this information is presently not known. The interpretation of the last form factor d\(t) in Eq.([T]) 
is less trivial but of equal significance for understanding the nucleon structure. It provides information on how the 
strong forces are distributed and stabilized in the nucleon [l7l . [l8j . 

In all theoretical approaches where it was studied so far, di = di{0) was found negative. This feature is expected to 
be deeply rooted in the spontaneous breakdown of chiral symmetry jl9U2ll. The form factor d\ (t) can be extracted 
from the beam charge asymmetry in deeply virtual Compton scatterin g 1201] . The EMT form factors of the nucleon 
have been studied in a variety of theoretical approaches: in lattice QCD I22|-[28l. in chiral perturbation theory [30l - [34j ]. 
in the chiral quark-soliton model [3oT - |42l | as well as in the Skyrme model|43|. Those of nuclei have also been studied 

It is also of great importance to understand how the nucleon undergoes changes in nuclear matter. Studying the 
EMT form factors of the nucleon in medium offers a new perspective on the internal structure of the nucleon, and an 
important step towards the understanding of how nucleon properties are modified in nuclei. The first experimental 
study of deeply virtual Compton scattering on (gaseous) nuclear targets (H, He, N, Nc, Kr, Xe) was reported in (47j . 
The admittedly sizable error bars of this first experiment did not allow to observe nuclear modifications. More precise 
future experiments of this type can in principle provide information on nuclear modifications of EMT form factors. 

Thus, in the present Letter, we aim at investigating the form factors of the total EMT of the nucleon in nuclear 
matter within the framework of an in-medium modified SU(2) Skyrme model, extending the previous work [43j . The 
Skyrme model [48|, |4|| provides a simple framework for the nucleon and connects chiral dynamics to the baryonic 
sector explicitly. Hence, the model can be easily extended to nuclear matter, modifications of chiral properties of 
the pion being taken into account. The changes of the nucleon in nuclear matter have been already examined within 
the in- medium modified Skyrme model [Hoi . Ifjlj . In Ref. (52j , the model has been further elaborated, the stabilizing 
term being refined in medium, which we will take as our framework to study the EMT form factors of the nucleon in 
nuclear matter. 
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2. Wc start with the in- medium modified chiral Lagrangian [5j 



£* = Tk THdoUdW) ~ a p ^r Tr(VC7) • (W f ) + — L- Tr^^tf), ^(cW)] 2 + a s ^|^ Tr (^ - 2 ) - ( 2 ) 



32e 2 7 



where {/ denotes the SU(2) pion field, F T the pion decay constant, e a dimensionless parameter, and m„- the pion 
mass. The medium modifications are contained in the following functions [50M5H ] 



a p {p) 
a s (p) 
70) 
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The a s ^ p depend on the S- and P-wave pion-nucleus scattering lengths, volumes (60 and en), and the density p 
of nuclear matter [50l [53j . and 5' is the Lorentz-Lorenz or correlation parameter. Near the threshold, the energy 
dependence of bo and Co being ignored and nuclear matter being approximated to be homogeneous, a StP simply 
become functions of the given nuclear density. Similarly, the function 7 in Eq. (JS|) describes the modification of the 
Skyrme term in nuclear matter proposed in Ref. [52j with 7 num and 7d C n fitted to the coefficient of the volume term in 
the semiempirical mass formula. We can treat the modified chiral Lagrangian in terms of the renormalized effective 



constants F* , 



F F* 



7 



1 / 2 e, and m* 



J 1 / 2 ™^. The behavior of these parameters 
in nuclear matter is consistent with those in chiral perturbation theory [54fl and QCD sum rules |55j . 

Homogeneous nuclear matter allows us to keep the hedgehog Ansatz for the pion field, i.e. U = exp [iT-F(r)] with 
a profile function F(r) in contrast to the case of local-density approximations for finite nuclei (5(| |57| . Consequently, 
one can immediately write the classical mass functional as 



M s * ol [F]=47r / dr 
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F 



sin 2 F /sin 2 P 



2e* 



2F' 



m *2p*2 



7T 7T.S 



(1 - cosF) 



(6) 



where F' = dF/dr. The minimization of the mass functional ([6]) leads to the equation for F(r) as follows 

„2 oip „;„2 



2Psin z F\ rF' F' 2 sm2F sin2F sin 2 Psin2F m; 2 r 2 sinF_ 

e* 2 Fl\ ) + ~ + e* 2 F:l 4 e* 2 Fl\r 2 4 "° 



(7) 



with the boundary conditions F(0) = tt and F(r) — > as r — > 00 imposed by the unit topological number of the chiral 
soliton. Having performed a collective quantization, we arrive at the modified collective Hamiltonian 
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where J 2 and I 2 are the squared collective spin and isospin operators, respectively, which act on the nucleon or A wave 
functions obtained from the diagonalization of the collective Hamiltonian (49|, and 8* is the moment of inertia of the 
soliton. A consistent description of the EMT form factors requires either to minimize the energy functional including 
rotational corrections, or to consider for the nucleon mass and other observables only the leading contribution in 
the limit of large number of colors N c [43[. In this work we will follow Ref. [43| and choose the second option. In 
particular, this means that in our treatment the nucleon and A masses (in vacuum or in medium) are degenerate and 
given by the minimum of the mass functional (JB]) . 

Input parameters in the Skyrmion sector are fixed as to t = 135 MeV, F v = 108.78 MeV, e = 4.854 following 
Ref. [IH . Those relevant to nuclear matter are determined by reproducing the coefficient of the volume term in the 
semiempirical mass formula and the experimental data for the compressibility of nuclear matter (52[: bo = —0.024 to" 1 , 
co = 0.09 to~ 3 , g 1 = 0.7, 7 nU m = 0.797 m~ 3 and 7d cn = 0.496 m~ 3 . All observables will be given as functions of p/po 
with normal nuclear matter density po = 0.5to 3 . 



From now on, the asterisks in expressions indicate the medium modified quantities which depend on the medium-dependent functions 
explicitly. Otherwise, we use the same symbol without any asterisk. 
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3. We are now in a position to calculate the EMT form factors of the nucleon in nuclear matter. Since the details of 
the general formalism have already been presented in Ref. [43j in free space, we briefly recapitulate only the necessary 
formulae here. The components of the static EMT [llj are given as follows: 
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(10) 



where T * is called the energy density. The density of the angular momentum pj, the pressure density, and the shear 
force density are expressed respectively as 
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The three form factors in Eq .([!]) arc finally given in the large N c limit as 
t 



(1 - cosF) 



Mm 



5M* N 2 



dm = 
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15M^ 
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d 3 r p*(r) 
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ry/-t 

where jo(z) and ji(z) represent the spherical Bcsscl functions of order and 1, respectively. At the zero momentum 
transfer t = 0, M|(0) and J*(0) are normalized as 

Mm = ^ |d 3 r T* (r) = 1 , J*(0) = | d 3 r p}(r) = 1 . (17) 

These relations were proven in [43j for a Skyrmion in free space, and we have checked that they hold also in nuclear 
matter. It is a feature of the approach with homogeneous nuclear matter that all expressions and proofs are formally 
analogous to the free nucleon case. In our approach, only the parameters of the Skyrme model are modified, but not 
the structure of the Lagrangian. (Of course, for a free nucleon the EMT is a Lorentz tensor. In nuclear medium 
we deal with a conserved Noether current.) Another important relation for p*(r) that is a consequence of the EMT 
conservation is the stability condition 



drr 2 p*(r) = . 



(18) 



Again we verified that the analytic proofs of Eq. (|18l) formulated in [43j for a Skyrmion in free space can be carried 
over straightforwardly to the medium situation. Finally the constant d\ in nuclear matter is given by 



d\ = 57r Af^r / dr r A p* (r) 



4ttM* n 



dr r 4 s*(r) 



(19) 



'o 3 j 

The conservation of the EMT implies that p*(r) and s*(r) arc connected to each other by the differential equation 
| s*'(r) + - s*(r) +p*'(r) = which is at the origin of the equivalent expressions for d* in (p~9|) . We have checked that 



also in the medium-modified Skyrme model this differential equation holds, and the different expressions in Eq. (|19j) 
yield the same result for d\, which is another demonstration of the theoretical consistency of the approach. Thus, the 
shear force density is directly related to the pressure density, so that we will concentrate only on the pressure here. 

In the chiral quark-soliton model [4l| it was observed that the negative sign of d\ is a consequence of stability, which 
requires p(r) > (repulsion) in the inner region, and p(r) < (attraction) in the outer region, which must balance 
each other according to Eq. (|T5)) . By inserting a factor of r 2 in the stability integral (fT8")) one basically obtains the 
expression for d\ in Eq. (fT9)) in terms of p(r). The "additional factor" r 2 gives more weight to the outer reg ion and 
is responsible for the negative sign of d\ . This pattern was also observed for the Skyrmion in free space [43| . Below 
we will recover this picture, with important medium modifications. 
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4. We now proceed to present the results of the EMT form factors of the nucleon in nuclear matter and discuss their 
physical implications. In Table |H we list the quantities relevant to the EMT form factors and their densities. First 
we want to mention that the energy density in the center of the nucleon is rather sensitive to the input parameters in 
the Skyrmion sector. In fact, the energy density in the center of the nucleon in free space in this work turns out to be 
T oo (0) = 1.45 GeV • fm~ 3 , while T oo (0) = 2.28 GeV • fm~ 3 was obtained in Ref. 0]. The reason lies in the different 
parameter set we use in the present approach. 

As the density of the surrounding nuclear medium increases from zero to normal nuclear matter density, the energy 
density in the center of the nucleon decreases by about a factor of 2. At the same time we observe that the mean 
square radius for the energy density 

Wff ,20) 

increases. This implies that the nucleon is swollen in nuclear matter and its energy density spreads to larger distance 
in comparison with that in free space. 

Let us more closely look into the change of the energy density in nuclear matter. The left panel of Fig. [T] shows 
how the energy density Too(j") changes in nuclear matter. We present in fact 4-7rr 2 T c J u (r) normalized by the effective 
in-medium nucleon mass such that the integration of the curve in the left panel of Fig. Q] yields unity. The solid 
curve draws it in free space, while the dashed and dotted ones depict those at p = 0.5 po and at normal density 
(p = po), respectively. Note that the nucleon mass also changes as the density does. As p increases, the energy 
density gets weaker and shifted to the large distance. Moreover, at a finite p, it falls off more slowly than that in free 
space, as the distance r gets larger. 



(rlo) 



TABLE I: Different quantities related to the nucleon EMT densities and their form factors: Tq (0) denotes the energy in 
the center of the nucleon; (tqq)* and {r 2 j)* depict the mean square radii for the energy and angular momentum densities, 
respectively; p*(0) represents the pressure in the center of the nucleon, whereas Tq designates the position where the pressure 
changes its sign; d* is the value of the dl(t) form factor at the zero momentum transfer. 
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FIG. 1: (Color online) In the left panel, the energy densities of the nucleon normalized by the nucleon mass, 47rr 2 Too (r)/Mj^, 
are presented as functions of r. The solid curve draws that in free space, whereas the dashed and dotted ones depict those 
at the density p — 0.5 po and at normal density po, respectively. In the right panel, the densities of the angular momentum 
normalized by the nucleon spin Sn = 1/2, 4-7rr 2 p*j(r)/ Sn are rendered as functions of r. The other notations are the same as 
in the left panel. 
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FIG. 2: (Color online) The pressure densities p*(r) and 47rr 2 p*(r) as functions of r in the left and right panels, respectively. 
Notations are the same as in Fig. [1] 



The density for the angular momentum p}(r), which is related to T * fe , vanishes in the center of the nuclcon. The 
corresponding mean square radius (rj)* which is denned analogously to (^qq)* in Eq. (|20|) starts to increase mildly 
as the density of nuclear matter increases as shown in Table U and exhibits similar behavior to (rg )*. 

Knowing the pressure density is also of great significance, which contains information on the spatial components of 
the EMT. In the fifth column of Table HI we list the values of the pressure in the center of the nucleon, p(0)* , given 
nuclear matter densities. It decreases with p increased. In the meanwhile, the position Tq, where the pressure changes 
the sign, becomes larger when p is switched on, which is in line with the features of the mean square radii discussed 
previously. That is, one can say that the shape of the nucleon is bulging out with nuclear matter. 

Figure [2] shows how the pressure is distributed in the nuclcon. As discussed already in Ref. [43| , the positive sign 
of the pressure for r < ro indicates repulsion, while the negative one in the region r > ro signifies attraction. In the 
Skyrmc model, the repulsive part is provided by the 4-derivative stabilizing (Skyrme) term, whereas the attractive 
one mainly comes from the kinetic term. Since the coefficient e is related to the vector meson coupling [59| . such a 
distribution of the pressure can be interpreted as follows: While its repulsive part (or core part) is mainly governed 
by the vector mesons {p meson), the attractive part (or long-range part) is explained solely by the pions. This picture 




r [fm] r [fm] 



FIG. 3: (Color online) The decomposition of the pressure densities 47rr 2 p*(r) as functions of r, in free space (p = 0) and at 
p — pa, in the left and right panels, respectively. The solid curves denote the total pressure densities, the dashed ones represent 
the contributions of the 2-derivative (kinetic) term, the long-dashed ones are those of the 4-derivative (stabilizing) term, and 
the dotted ones stand for those of the pion mass term. 
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FIG. 4: (Color online) The EMT form factors of the nucleon M% (t), J*(t), and d*(t) as functions of t for the nucleon. Notations 
are the same as in Fig. [T] 



is right also in nuclear matter, which is explained in Fig. [3] Comparing each contribution in the right panel (p = po) 
to the corresponding one in the left panel, we immediately realize that the 2- and 4-dcrivative terms are suppressed 
as well as stretched to larger distance. On the other hand, the mass term is enhanced at finite density. In order to 
understand this feature more precisely, we examine the expression for p*(r) given in Eq. (|12[) . The relevant coefficient 
in the kinetic term is F '* s that becomes smaller as p increases, so that its contribution to p*(r) gets diminished. In 
the meanwhile, e* governs the strength of the contribution of the stabilizing term. It grows as p is turned on, which 
leads to lessen the contribution of the stabilizing one to p*(r). In the case of the pion mass term contribution to 
p*(r), the pertinent coefficient is m* 2 F* 2 s = m 2 F 2 a s , which becomes larger with p increased. That results in the 
enhancement of the pion mass term. While all these contributions undergo changes in nuclear matter, the stability 
condition of Eq. (|T5|) is still satisfied. Integrating each contribution, we obtain for the case of free space: 



f 

Jo 



-9.550 McV for the kinetic term 
drr 2 p(r) = { 12.446 MeV for the Skyrme term (21) 
-2.896 McV for the pion mass term , 



whereas for the case of p = po 



f 

Jo 



-7.741 MeV for the kinetic term 
drr 2 p*(r) = { 13.652 McV for the Skyrme term (22) 
—5.911 MeV for the pion mass term . 



We see that adding up all contributions becomes zero, as required by the stability condition (fl~8"|) . 

Having performed the Fourier transforms of the densities T * (r), p*j(r), and p*(r) discussed so far, we immediately 
obtain the nucleon EMT form factors. Figure [4] draws the resulting EMT form factors as functions of t. The form 
factors M|(£) and J*(t) fall off more rapidly as p increases. On the other hand, d\{t) is rather distinguished from 
the other two form factors. While M|(i) and J*(t) are constrained to be 1 and 1/2 at t = as shown in Eq. (|17|) . 
d\{t) does not have such a constraint. However, had one scaled d\(t) so that their values may coincide at t = 0, d\(t) 
would have ended up with almost the same response to p as M% (t) and J* (t) . As shown in Fig. HI all the EMT form 
factors get stiffer near the zero momentum transfer, as the density of nuclear matter increases. It implies that all the 
corresponding radii get enlarged with p increased, as was already discussed for the case of (t~ 2 )* and (r 2 j)*. 

All EMT form factors can be well approximated by dipolc-type parametrizations F*(t) = F*(0)/(l — t/Mp 2 ,) 2 with 
the generic dipole masses Mp» listed in Table HH 



TABLE II: The dipole masses of the EMT form factors M| (t), J*(t) and d*(i) which approximate well the results in Fig. [4] 
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5. So far we have considered the medium effects in the range < p < po where our approach based on linear 
response theory of pions in medium [53j is applicable. The description in terms of hadronic degrees of freedom is 
expected to become inappropriate and to break down at densities beyond pq. It is interesting to remark that by 
continuing our in-medium modified Skyrme model to p > po we find indications for this expected breakdown. We 
would like to stress that we do not expect the naive extrapolation of the linear response in Eqs. (|3j H]) to work 
quantitatively. However, the results are qualitatively insightful. 
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FIG. 5: (Color online) In the left panel, the correlated change of p*(0) and Tq O (0) are drawn with p varied. In the right panel, 
the T(* /Mn and p depicted as a function of p/po- The maximal density is given as about 6.74 po, above which the Skyrmion 
does not exist anymore. The filled circle on the solid curve represents the value of Tqq/M^ at normal nuclear matter density. 



The left panel of Fig. [5] shows the correlated modifications of the pressure and energy densities in the center of the 
nucleon. We already have seen from Figs. Q] and [5J both p*(0) and T o * o (0) decrease as p grows. This trend continues 
also at higher densities, and the interesting point is that p*(0) and Tq O (0) are diminished at the same rate. This is 
natural if we consider nuclear matter density as an external parameter. The response of a stable system to variations 
of external parameters must be such that the pressure increases if the energy density does. The relation of p*(0) and 
Tq O (0) remains linear to a good approximation, until the Skyrmion ceases to exist at the "maximal nuclear density" 
Pmax = 6.74/50. This maximal value arises because a p (p) — > for p — » p m „*. In this limit the energy functional ([6]) has 
no minimum. This can be shown by means of the Derrick theorem [58j , which is equivalent to the stability condition 
(HU) 53. In other words, in the limit 

P Pmax the stability condition breaks down. We stress that this maximal 
density appears only in the mathematical description of our medium-modified Skyrmion, and has nothing to do with 
the physical critical point relevant to the restoration of chiral symmetry. 

The fact that the Skyrmion ceases to exist as p — > p max is of interest from a mathematical point of view. Of far 
greater relevance is that physical description of the nucleon in medium breaks down before this point is reached. This is 
because at p « (2-3)po two things happen. Firstly, the normalized energy density T o * o (0)/M^ (alternatively one could 
also consider the baryon number density) becomes comparable with the density of the surrounding nuclear medium 
(see the right panel of Fig. [5J . Secondly, the size of the nucleon, as measured for instance in terms of the square root 
of the mean square radii or rg, becomes comparable to p -1 / 3 which implies a spatial overlap of the nucleons in the 
nuclear medium. At this point it makes no sense anymore to mention an "individual nucleon" in nuclear medium, 
and a description in terms of more appropriate degrees of freedom becomes necessary ( "quark matter" ) . It is very 
interesting to observe that the medium-modified Skyrme model signals in this way its limitations. 



6. In the present work, we investigated the energy-momentum tensor form factors of the nucleon in nuclear matter, 
based on the in-medium modified Skyrme model. Employing all the parameters fixed in previous works, we derived 
the densities relevant to the energy, the angular-momentum, the pressure, and the shear-force distributions in the 
nucleon. We first examined how the energy and angular-momentum densities were modified in nuclear matter. In 
general their values in the center of the nucleon are suppressed but they are stretched to larger distance in nuclear 
matter. This was documented by the increase of the corresponding mean square radii with nuclear matter density. 

We also analyzed the medium modifications of the pressure density. Again, its value in the center of the nucleon 
decreases as nuclear matter density grows. The change of the pressure density showed a similar behavior to the energy 
density in matter. The position where the pressure changes sign increases with nuclear matter density increased, and 
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the pressure density falls off in medium more slowly than that in free space. Since the pressure plays an essential role 
in the stability of the Skyrmion, we scrutinized each contribution of the kinetic, stabilizing, and mass terms. With 
nuclear matter density increased, the contributions of the kinetic and stabilizing terms are suppressed in magnitude, 
while that of the mass term becomes larger, and all contributions extend to larger distances. The stability condition 
was shown to be satisfied also in nuclear matter. Moreover in nuclear matter the constant d\ is negative as a 
consequence of stability. In particular, we observe that medium effects further decrease the value of d\. 

The energy- momentum tensor form factors M2(t) and J(<), which are constrained to be 1 and 1/2 at zero momentum 
transfer, were found to fall off more rapidly as nuclear matter density increased. The absolute magnitude of the values 
of the d[(t) form factor increased in medium, compared to that in free space. The medium effects on its ^-dependence 
were more or less similar to the other two form factors. Our results may pave a way towards a better understanding 
of the discrepancy between predictions obtained from different models of nuclei [TJ H3, [45| . 

The present results for the modification of the EMT form factors and related properties, with nuclear density varied 
from zero to normal nuclear matter density and beyond that, indicate a breakdown of the description of the nucleon 
in nuclear matter densities at p (2-3)po- At such high densities the nucleon in nuclear matter starts to loose its 
identity, and it would finally melt away, which implies that a description in different degrees of freedom becomes 
appropriate at such high densities. 
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